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Qubits are neither fermions nor bosons. A Fock space description of qubits leads to a mapping 
from qubits to parafermions: particles with a hybrid boson-fermion quantum statistics. We study 
this mapping in detail, and use it to provide a classification of the algebras of operators acting on 
qubits. These algebras in turn classify the universality of different classes of physically relevant qubit- 
qubit interaction Hamiltonians. The mapping is further used to elucidate the connections between 
qubits, bosons, and fermions. These connections allow us to share universality results between the 
different particle types. Finally, we use the mapping to study the quantum computational power 
- - , of certain anisotropic exchange Hamiltonians. In particular, we prove that the XY model with 

(N . 

nearest- neighbor interactions only is not computationally universal. We also generalize previous 
\ results about universal quantum computation with encoded qubits to codes with higher rates. 

(N : 

I. INTRODUCTION 

. It is an experimental fact that there are only two types oi fundamental particles in nature: bosons and fermions. 
QQ ■ Bosons are particles whose wavefunction is unchanged under permutation of two identical particles. The wavefunction 
^SJ : of fermions is multiplied by —1 under the same operation. An equivalent statement is that bosons transform according 
to the 1-dimensional, symmetric, irreducible representation (irrep) of the permutation group, while fermions belong 
to the 1-dimensional antisymmetric irrep. The permutation group has only these two 1-dimensional irreps. What 
about particles transforming according to higher-dimensional irreps of the symmetric group? Much research went into 
studying thispossibility, in the early days of the quark model, before the concept of "colored" quarks gained widespread 

; acceptance JlJ,^]. However, there are now good reasons to believe that particles obeying such "parastatistics" do not 
exist (Ref. ^ p. 137]). Nevertheless, as we will show below, the traditional definition of a Hilbert space of qubits is 

, inconsistent with the properties of either bosons or fermions. 

' The description of bosons and fermions in terms of their properties under particle permutations uses the language 
\ of first-quantization. A useful alternative description is the second-quantized formalism of Fock space A basis 

■ state in the boson or fermion Hilbert-Fock space can be written as \n",n2, ■■■), where nf counts how many bosons 

Q_i' (a = b) or fermions (a = /) occupy a given mode, or site i. Note that the total number of modes does not need to be 
^li specified in the Fock-basis. Ignoring normalization, raising, aj (lowering, a^) operators increase (decrease) nf by 1. 

Ch a consequence of the permutation properties of bosons and fermions is that their corresponding raising and lowering 

^ ' operators satisfy commutation and anti-commutation relations: 

[6j,6]]=0, [6i,5]]=% bosons 
{//,/]} = 0, {h, /]} = % fermions. 

From this follow a number of well-known facts |^,^. Let nf — ajai; this is the number operator, which is diagonal in 
the Fock-basis \nf ,n2 , ■■■) , and has eigenvalues nf. Then: 

• [&J, fej] = =^ an arbitrary number of bosons nf can occupy a given mode i. On the other hand, {//, /|} = 



only n^ ~ 0,1 is possible for fermions. 

• [fe,j,fet] ^ g.^ ^ -t-i^g Hilbert space of bosons has a natural tensor product structure, i.e., \n\,n2, ...) = \n\) (g) 
1 712) (8) • • •. More specifically, it is possible to independently operate on each factor of the Hilbert space. However, 

{/.,/]} = S.,, =^ ...,n^^_i, 1,4+1, ...) = {^l)J:r>i\n(, ...,n^_„0,nf^„ ...), 

which means that the outcome of operating on a mode of a multi-fermion state depends on all previous modes 
(the order of modes is actually arbitrary). This non-local property means that the fermionic Fock space does not 
have a natural tensor product structure, although it can be mapped onto one that does using the Jordan- Wigner 
transformation [|j (see ||] for a more detailed discussion). 

What about qubits? The standard notion of what a qubit is, is the following |^: 
Qubit: 
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• A qubit is a vector in a two-dimensional Hilbert space Hi = span{|0)i, (like a fermion), and 



• An A^-qubit Hilbert space has a tensor product structure: Ti. — '^iLiHi (like bosons). 

It appears that a qubit is a hybrid fermion-boson particle! We conclude that qubits do not exist as fundamental 
particles. This motivates us to consider an intermediate statistics of "parafermions" in order to have a Fock space 
description of a qubit. We define the parafermionic commutation relations by 

{a^,al} = 1, 

[a„a]]-0 if i^j. (1) 

Here i,j are different modes, or different qubits. The relation [ai,a]] = for i ^ j immediately implies a tensor 
product structure, while {ai,a|} = 1, which together with ai\0) — (|0) is the vacuum state) implies 

aitti = alal = (2) 

in the standard (irreducible) two-dimensional representation. Therefore a double-occupation state cannot be realized, 
i.e., the single-particle Hilbert space is two-dimensional. These are exactly the requirements for a qubit. 

In fact, the notion of particles with "intermediate" statistics such as parafermions is well known and established 
in condensed matter physics, e.g., hard-core bosons, excitons, or the Cooper pairs of superconductivity pCj] (see also 



Section VI). Such particles are always composite, i.e., they are not fundamental. Another way of obtaining a particle 
that is neither a boson nor a fermion is to simply ignore one or more degrees of freedom. This is by and large the 
approach taken in current proposals for the physical implementation of quantum computers. For example, a single 
spin-i, without the orbital component of its wavefunction, behaves exactly like a qubit. This is the case of the 
electron-spin qubit in quantum dots |pT| . Related to this, a truncated multi-level atom can also approximate a qubit, 
as in the ion-trap proposal What are the implications of this for quantum computing (QC)? In a nutshell, "ideal" 
qubits are hard to come by. If a qubit is to exist as an approximate two-level system, or as a composite particle, 
or as a partial description of an object with additional degrees of freedom, this means that some robustness is lost 
and the door is opened to decoherence. E.g., the additional levels in a multi-level Hilbert space can cause "leakage", 
the orbital degrees of freedom act as a bath coupled to the spin-qubit, and a composite particle may decay (e.g., the 
exciton-qubit [p^). 

The advantages of the parafermionic formalism for qubits, however, are not necessarily in understanding these 
sources of decoherence, because this formalism "accepts" qubits as particles. Instead, the parafermionic formalism 
allows us to naturally establish mappings between qubits, fermions, and bosons. This mapping serves to transport 
well-known results about one type of particle to another, which, as we show below, clarifies questions regarding the 
ability of sets of one type of particle to act as universal simulators jlj] of sets of another type of particle. It also helps 
in connecting the Hamiltonians of condensed matter physics to standard tools of quantum computation. 

The structure of the paper is as follows. In the next section we formally introduce the second quantization of qubits. 



We then classify the algebraic structure of parafemionic operators in Section [II . This classification, into subalgebras 
with different conservation properties, is very useful for establishing which subsets of qubit operators are universal, 
either on the full Hilbert space, or only on a subspace. This is taken up in the next two sections, where we establish 
the connection between parafermions and fermions (Section ^) and bosons (Section ^). The connection to fermions 
and bosons also works in the opposite direction: we are able to classify which fermionic and bosonic operator sets are 



universal. This has implications, e.g., for the linear optics quantum computing proposals |13,n6|. Section VI shows 



how to construct parafermions out of paired fermions and bosons, emphasizing the compound-particle aspect of qubits 



With the connections between fermions, bosons, and parafermions clarified, we explain in Section VII a remarkable 



difference between parafermions and the other particle types: bilinear parafermionic Hamiltonians are sufficient for 



universal quantum computation, whereas fermionic and bosonic Hamiltonians are not. In Section VIII we briefly use 



the mapping to fermions to derive the thermal fluctuations of non-interacting parafermions at finite temperature. In 



Section IX we apply the classification of the various parafermionic operator subalgebras to the problem of establishing 
universality of typical Hamiltonians encountered in solid state physics. We generalize a number of our previous results 
In particular, we establish that the XY model is not universal with nearest-neighbor interactions only; and, 
we prove universality of the XXZ model for codes with arbitrarily high rates. We conclude in Section 

II. SECOND QUANTIZATION OF QUBITS 

As in the cases of bosons and fermions, a parafermion number operator in mode i can be defined as 
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Hi = alai, 

with eigenvalues Ui = 0, 1. The total number operator is n = X^i^*- ^ normalized basis state in the parafermionic 
Fock space is 

\---n^---)=l[ialr\0), 

i 

which we think of as representing a state with the i^^ qubit in the "up" ( "down" ) state if the i^^ parafermion is present 
(absent), i.e., = 1 (0). Qubit computational basis states are thus mapped to parafermionic Fock states. Equivalently, 
consider the following mapping from qubits to parafermions: 

|Oi---0,_iO,0,+i---) ^ |0) 
|Oi---0,_ilA+i---> ->al|0), 

where on the left and 1 represent the standard (first-quantized) logical states of a qubit. Qubits can thus be identified 
with parafermionic operators. 

The mapping of qubits to parafermions is completed by mapping the Pauli matrices u" to parafermionic operators: 

CT+ a] a~ Oi erf ^ 2ni- 1. (3) 
It is then straightforward to check that the standard sl{2) commutation relations of the Pauli matrices, 

[a+,a-]=2%a,r 

are preserved, so that we have a faithful second quantized representation of the qubit system Hilbert space and 
algebra (Of course we could also have mapped su{2) = {cr^, cr^, cr^} to the parafermionic operators, by appropriate 
linear combinations.) To illustrate the multi-qubit Hilbert-Fock space representation, consider the case of two modes, 
i.e., i,j = 1,2. The space splits into a vacuum state |00) = |0), single-particle states |01) = a\ |0) and |10) = a^ and 
a two-particle state |11) — a\a2 |0). It is important to emphasize that the parafermionic formalism is mathematically 
equivalent to the standard Pauli matrix formalism. We will be using both in the sections below, starting with the 
parafermionic, as it makes particularly transparent the translation of known results about fermions to qubits. 



III. GENERAL PROPERTIES OF PARAFERMIONIC OPERATORS 



A^-qubit operators in QC are elements of the group U{2^). We will begin our discussion by identifying a set of 
infinitesimal parafermionic generators for U{2^). Recall that with any r-parameter Lie group there are associated r 



infinitesimal generators |19|. E.g., in the case of su{2) these are, in the two-dimensional irreducible representation, 
the Pauli matrices {(^XT^^yT'^z}- Now, let a = {ai},/3 = where ai,j3j can be or 1. In terms of parafermionic 

operations, any element of U{2^) can be written as U{b) = eyip{~'iJ2a i3^'^^Qa.i3{N))i where b°'^ are continuous 
parameters (generalized Euler angles) and the 2^ x 2^ infinitesimal group generators Qa.pi^) Stie defined as follows: 

qUN^)^{a^j,r- ■■■{a\r\ qp{N - N^) ^ a^^^^ ■ ■ ■ a^,\ (4) 



let Na = Yh=i and 



Then: 

Qo^AN) - qi{N^)qp{N - N^). (5) 

The Qa.p{N) will be recognized as all possible transformations between A^-qubit computational basis states, e.g., for 
N = 2 the set of 16 operators is: 

{/, a|, Cj, ai, a2, a\a\, 0102, a|ai, a|a2, a\ai, a\a2, a\a\ai, a\a\a2, a\aia2, aja2ai, a\a\a2ai \ , where / is the identity op- 
erator. The set Qa,oiN) generates all possible basis states from the vacuum state. Hermitian forms are Q + and 



i{Q — Q^). We will turn to the hermitian set of generators in the discussion of applications, in Section IX. 

Note that infinitesimal generators are not the generators one usually considers in QC. Rather, in QC, a gate 
operation is obtained by the unitary evolution generated through the turning on/off of a set of physically available 
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Hamiltonians {H^}, that are generally a small subset of the 2^ x 2^ infinitesimal generators Qa,p{N). "Generated" 
here has the usual meaning of allowing linear combinations and commutation of Hamiltonians. We will say that a 
set of Hamiltonians is universal with respect to a Lie group Q if it generates the Lie algebra of that group. The 
question of the dimension of the universal set of Hamiltonians with respect to U{2^) is somewhat subtle, since it 
is context dependent. Lloyd showed that given two non-commuting operators A,B, represented by n x n matrices, 
one can almost always generate U{n) j2^. However, it is not necessarily clear how this result is related to physically 
available Hamiltonians, since in practice one may have only limited control over terms in a Hamiltonian. E.g., the 
standard Hamiltonian generators for SU{A) (two qubits) is the 5-element set {u^, (t|, cr^, erf, (tJ^(t|}. However, the 
4 -element set {erf , erf, erf erf — crfcr|,(Ti • 0^2} also generates 5'f/(4), and may be physically available iQ. Another 
example are the following sets of, respectively, five, four, and three generators: {erf , erf , erf , erf, erf erf }, {erf, erf , cierf -|- 
C2cr2 1 o'i'^'IIj Wi^ f 2 J cierf -I- C2erf -I- Caerf erf } (where c, are constants). Which set of generators is physically available 
(i.e., directly controllable) depends on the specific system used to implement the quantum computer. As we will show 
below, it is sometimes the case that a given, physically available, set of Hamiltonians is universal with respect to a 
subgroup of U{2^), which may be quite useful, provided the subgroup is sufficiently large (typically, still exponential in 
TV). This notion of universality with respect to a subgroup is what gives rise to the idea of encoded universality [pl|-p4|: 
one encodes a logical qubit into two or more physical qubits, and studies the universality of the subgroup-generating 
Hamiltonians with respect to these encoded/logical qubits. 

The infinitesimal parafermionic generator Qa.p{N) can be rearranged into certain subsets of operators with clear 
physical meaning, which we now detail. 

1. Local subalgebras: The tensor product structure of qubits is naturally enforced by [ai,a^|] = for i 7^ j. This 

induces a tensor product structure 0,^^ sli{2) on the subalgebras formed by the grouping sli{2) = {ai, a|, 1 — 
2ni}. Each sli{2) can only change states within the same mode. 

2. SAp - Subalgebra with conserved parity: Define a parity operator as 

It has eigenvalues 1 (—1) for even (odd) total particle number. The operators that commute with the parity 
operator form a subalgebra, which we denote by SAp. Let k (l) be the number of aj (a^) factors in Qa,i3{N), 
i.e., 

SAp consists of those operators having k — I even, so its dimension (i.e., number of generators) is 2^^/2. To see 
this, let Qi be in SAp, and consider its action on a state with an even number of particles \n). Since k — I is 
even, Qi\n) — \n') where n' is also even. Now, pQj\n) — p\n') = +\n'), but also Qip\n) = Qi{+\n)) = \n') so 
[Pj Q] = 0, e.g., for = 2 SAp consists of: {/, ala\, aia2, ajai, a\a2, alai, ala2, a2a\a2ai} . 

3. SAn - Subalgebra with conserved particle number. This subalgebra, which we denote SAn, is formed by all 
operators commuting with the number operator n. These are the operators for which k = I, so its dimension 

is SfeLo (^) = ^n\n\ • ^'^iS' ^® SAn, and consider its action on a state |n) with n particles. 

Qii cannot change this number since k = I, but it can transform \n): hQii\n) — n\n)' — n\n)' . However, 
Qiin\n) = nQii\n) = n\n)' , so [Qu, h] — 0, e.g., for N — 2 SAn consists of: {/, a\ai, a\a2, a\ai, a\a2, a\a\a2ai \ . 
Clearly, SAn cSAp. 

4. Subsets of bilinear operators: There are two types of bilinear operators for i ^ j: ajaj (which conserve the 
particle number), and aiaj,a|aj (which conserve parity). Let /i = {ij), then first: 

= a] a, + a] a J 

= n., - n, (6) 
and Ty = i[TI^,T^] form an su{2) subalgebra, denoted su^(2). Clearly, sit*,(2) eSAn. Second: 

— QiUj + a J a] 

Rl^n,+ rij - 1 (7) 
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and form another su{2) subalgebra, denoted su^(2) sSAp. Note that [su^(2), suj^(2)] = since any product 

of raising/lowering operators from these algebras contains a factor of a^ai or a\a\. Consider as an example the 
case of iV = 2 modes. Whereas the direct product group SUi{2) ® SU2{2) yields all product states, the group 
SU*{2) © SU'^{2) can transform between states with equal particle number and states differing by two particle 
numbers. 

5. Generators of SAn(7V): The set of Hamiltonians {ajaj}f t\ generates SArt(7V), i.e., the subalgebra of conserved 



particle number on N modes (qubits). Proof: this set maps to the XY model (see Section [IX B| ). The rest follows 
using the method of [Q. Note that {ajaj}fj=i does not generate SAn(A^ + 1), since this set cannot generate 
hih2 ■ ■ ■ fiN- 

6. Generators of SAp{N): The set of Hamiltonians {ajoj, UiUj + aja], i{aiaj — alaj)}fj^i yield all states with even 
particle number on N modes from the vacuum state. (Proof is trivial.) 

7. Generators of SU{2^): In order to transform between states differing by an odd number of particles it is 
necessary to include the operators {a^, a|} as well. The corresponding set {ajaj, OiOj, a|aj, a^, al}^j^^ generates 



a set of universal gates (proof is trivial) , and then by standard univer sality results [g5U26| the entire SU{2^) 



Additional structure emerges from a mapping between fermions and parafermions. This structure can help us both 
in simulating fermionic system using qubits, and the understanding of universality of qubit systems. 



IV. FERMIONS AND PARAFERMIONS 

A general fermionic Pock state is 

\ni,n2,- ■ ■)f ^ (8) 

where = 0, 1 is the occupation number of mode i. As is well known p7| , the fermionic ("supergroup" |]l^ ) U(2^) 
has infinitesimal generators 

w = (/^r" • • • (/Ir^/^" • • • /f^ 

where 

N 

A = (g)(l-n.)- 

This basis is equivalent by a linear transformation to the more familiar set 

QiA^) = ifhr'' ■■■{fir fN 

which transforms between all possible fermionic Pock states ( "fermionic computational basis state" ) . There is a group 
chain of this group, 

U{2^) D S0{2N + 1) D SO{2N) D U{N) (9) 



and the generators of the subgroups are known |19 



The Jordan- Wigner ( JW) transformation Q , recently generalized in Ref . |2^ , allows one to establish an isomorphism 
between fermions and parafermions. Defining 

S/^(g)(l-24), 5,^(g)(l-2nfc), (10) 

k=l k=l 



the mapping is: 



n{ -> m, 

fi ^ 0-iSi, 

fl - 4S^. (11) 
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The action of the fermionic operators on the state (g) is equivalent to that of the corresponding parafermionic operators 
on the state |ni, n2, • • •). To see this, note that [a^, Si] = 0. Therefore the effect of the JW transformation is quite 
simple: by commuting all Si to the left when when mapping a fermionic infinitesimal generator to a parafermionic one, 
we see that (i) the parafermionic a^, a] operators will yield a state with the same parafermionic occupation numbers as 
the corresponding fermionic state, (ii) the action of the product of S'^'s is to produce a phase ±1. (This may become 
a relative phase when acting on a state that is a superposition of computational basis states.) This allows us to study 
algebraic properties of one set of particles in terms of the other. 

Using the JW transformation we find that the same subgroup chain holds for parafermions, and we can im- 
mediately write down also the infinitesimal generators for the corresponding parafermionic subgroups. The result is 
given in Table ^ 

The significance of these subgroups for QC is in the classification of the universality properties of fermionic and 
parafermionic Hamiltonians. E.g., a Hamiltonian of non-interacting fermions, i.e., one including only bilinear terms 
{fi fj: fifjT fj fi} is not by itself universal since it merely generates SO{2N). Recent work has clarified what needs 
to be added to such a Hamiltonian in order to establish universality ||6|, p9| , ^ . Regarding S0{2N + 1), note that one 
must carefully discuss the hcrmitian terms ft + fj and i{fi — fj) if one wants to consider them as Hamiltonians, since 
it is unclear which physical process can be described by such Hamiltonians (a single fermion creation/annihilation 
operator can turn an isolated fermion into a boson, a process that does not seem to occur in nature). 

A more powerful classification, from the QC viewpoint, is in terms of physically available Hamiltonian generators 
of the subgroups. An interesting restriction of the set of infinitesimal generators to a physically reasonable set of 
Hamiltonians is to consider only nearest-neighbor interactions, where possible. The results known to us in this case 
are presented in Table ||. 

A couple of comments are in order regarding Table ||: First, note the group SO(2N + 1) may be unphysical not 
just for fermions since its generators must contain terms like fi + fj in its Hamiltonian, but also for parafermions: it 
requires a non-local Hamiltonian due to the Si term. Second, the corresponding fermionic generators for U{2^) given 
here is unphysical because it includes terms that are linear in fi and furthermore non-local. A physically acceptable 
set is {fj fi+i, fifi+i, fi fi^ififi+i, h.c.}, but this set is not universal over the full 2^-dimensional Hilbert space 
(since it conserves parity). This means that a qubit needs to be encoded into two fermions in this case, a situation 
we explore further in Section Now let us verify the claims of Table Our strategy is to show that in each case, 
we can use the Hamiltonians for generating all infinitesimal generators of the corresponding subgroup in Table |. 

Consider first the subgroup SU{N): In the fermionic case, we claim that this subgroup has nearest neighbor 
Hamiltonian generators fjfi+i and their hermitian conjugates. E.g., for iV = 3, if we have the four operators /1/2, 
flfs and h.c, then we can generate flfs = [/I/2, flfs] and h.c, as well as n{ - = [///j , /]/*]• This yields a total 
of nine operators, eight of which are linearly independent, that generate SU{3). As for parafermions, we can use the 
JW transformation to get fj^ifi al_^_-f^Si+iaiSi = a|^i(l — 2hi)ai = al_^_^ai (where we have used [0^,6'^] = and 
fiiQi = a\aiai = 0). This establishes an isomorphism between the fermionic and parafermionic generators for SU{N). 
Hence the parafermionic subgroup SU{N) is generated by a]^aj+i and h.c. 

Now consider SO{2N): In the fermionic case we have /|/|, and using the result for U{N) we also have /j/i; 
therefore we have [/l/ij/^/J^] = /l/l -Clearly, the interaction range can be extended to cover all generators. For the 
parafermionic case, using the JW transformation we find fj^ifj al_^-^^Si+ialSi = a|^^i(l — 2ni)al = al_^^al, so that 
we again have an isomorphism with the fermionic case. 

Next consider the (unphysical) subgroup SO{2N+l): In the fermionic case it suffices to note that fj] = fifj and 

^{fi' fj] ~ fi fj' that we can generate all infinitesimal generators by the linear terms fi and fj . The parafermionic 
case follows by the JW-transformation. 

Finally, in the U{2^) case the universality of the parafermionic set {a^, a\ai+i, h.c.} follows from that of the set of 
all single qubit operations together with the Hamiltonian of the nearest-neighbor XY model [Eq. (|l^) below] , proved 
in Ref. ||3^. The fermionic case follows by the JW-transformation. 

Let us recapitulate the meaning of the results presented in this section: we have shown how to classify subalgebras 
of fermionic/parafermionic operators in terms of the groups they generate. This therefore classifies their universality 
properties with respect to these groups. This is particularly important in the context of a given set of physically 
available Hamiltonians. Our method employed a mapping between fermions and parafermions, which allowed us to 
easily transport known results about one type of particle to the other. 
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V. BOSONS FROM PARAFERMIONS 



A linear combination of different-mode parafermions can approximately form a boson. Define 



f ^ 



B = —= > flj. 



Then using Eq. we have 

TV 
i=l 

If the paraferniion number is much smaller than the available number of sites/modes, i.e., when n <ti N, then 
[B, B^ w 1, which is an approximate single-mode boson commutation relation. 

To get K boson modes, we can divide TV into K approximately equal parts. Each part has Na = N/K qubits and 
approximately represents a boson. The k^^ boson is Ba — -^^= ■ Then 

[Bo,,B'p\^ 5ap[l- —) w dap. 

Physically, this means that a low-energy qubit system (with most qubits in their ground state) can macroscop- 
ically behave like a boson, or a collection of bosons. If the Hamiltonian is of the bilinear form H = —B^B 



~jj + O'iO'j^ , the ground state with n <ti N parafermions is (B^)^ |0), i.e., n [(-B^)" |0)] w n [(-B^)" |0)] . 

A separate conclusion that follows from this result is that a low-energy non-interacting qubit system can naturally 
simulate the dynamics of bosons. 



VI. PARAFERMIONS FROM FERMIONS AND BOSONS 



As stated in the Introduction, qubits do not exist as fundamental particles. This means that they are either 
approximate descriptions (e.g., a spin in the absence of its spatial degrees of freedom), or have to be prepared by 
appropriately combining bosons or fermions. I.e., a qubit can be encoded in terms of bosons or fermions under certain 
conditions (see also ||33)). We consider bosonic or fermionic systems with 2N single-particle states. Let fc = 1, 2, . . . , 
denote all relevant quantum numbers (including spin, if necessary). The following three cases yield parafermions. 

Case 1: Fermionic particle-particle pairs — Under the condition n2i^_^ = n^^. it can be shown that 
{/2fe/2fc-i, flk-ifL} = 1 ^^'^ [f2k-ihk, fli-jli] = for fc 7^ L Furthermore, the set {f2k-i, flk-ifL^ '^L-i +"-2fe ^ 1} 
satisfies the commutation relations of sl{2). Therefore the mapping <J=^ /2fc/2fc-i, o^. flk-iflk ^'^^ 

2nfc <;=4> fi2fc-i + "2fc' ^ mapping to parafermions. The vacuum state of parafermions in this case corresponds 
to the vacuum state |0)/ of fermions. Example: Cooper pairs. 

f f . 

Case 2: Fermionic particle-hole pairs — Under the condition tIj^.i + n2^, = 1 it can be shown as in Case 1 that 
Ok ■^^=> flkf^k-i, a[. -<=4> flf._if2k and 2nk — 1 <J=^ ^2fc-i ~ "-2/0 ^ mapping to parafermions. However, in this 
case the vacuum state of parafermions is |0) = f^j^ ■ ■ ■ f\f\ |0)^, because then |0) — for all k. This vacuum state 
plays the role of a Fermi level. Example: excitons. In fact, all quantum computer proposals that use electrons, e.g., 
quantum dots electrons on Helium jsj,^^, are equivalent to this case. E.g., /2/1 and f\f2 can represent the 
transition operators between two spin states in the quantum dot proposal. 

Case 3: Bosonic 'particle-hole' pairs — Under the condition n\j^_-^ + n\j^ = 1 it can be shown as in Case 1 that 

Ofc <J=^ b\jj]2k-i, flfe *^=^ ^2fc-i^2fc and 2nk — 1 "-2fc-i ~ "2/0 ^ mapping to parafermions. However, in this case 
the vacuum state of parafermions is |0) = ' ' ' ^\k ' ' ' ^4^2 |0)i.i again because then |0) = for all k. Example: 
dual-rail photons in the optical quantum computer proposal 

This classification illustrates the by- necessity compound nature of a qubit, and puts into a unified context the 
many different proposals for constructing qubits in physical systems. Note that it is possible to use more than two 
fermions or bosons to construct a parafermion. Further implications, especially as related to the simulation of models 
of superconductivity (Case 1) on a quantum computer, have been explored in Ref. [p6[. 
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VII. PARAFERMIONIC BILINEAR HAMILTONIANS ARE UNIVERSAL BUT FERMIONIC AND 

BOSONIC ARE NOT 



In this section we discuss an apparently striking difference between tlie universality of bilinear Hamiltonians acting 
on fermions and bosons, as compared to parafcrmions. Let us consider the set of particle-number-conserving bilinear 
operators of bosons, fermions and parafermions: 

bib J, f!fj, alaj. 

As noted in Table |, in the fcrmionic case these operators generate the group U {N) where N is the number of particles. 
The same is true for bosons . Clearly, therefore, fermionic and bosonic Hamiltonians containing only these operators 
are not universal with respect to an interesting (i.e., exponentially large) SU{2'^) subgroup. On the other hand, as 
discussed in the previous section, these fermionic and bosonic operators can be used to define parafermionic operators 
ajflj in two-to-one correspondence. As mentioned in Section the set {alaj}fjt\ generates the subalgebra SAn(iV), 

with dimension ^j^j (> 2^) (recall that the total number of Qa.piN) operators is 2^^). The corresponding Lie group 
appears to be large enough to be interesting for universal quantum computation. This expectation is borne out, 
since one can construct an XY model, Eq. (^ below, using the set {a\aj}. As shown in Ref. ||2^, the XY model 
is by itself universal provided one uses three physica l qubit s per encoded qutrit, together with nearest-neighbor and 



next-nearest-neighbor interactions (see also Section |IX D l| ) . We discuss the XY model in detail in Section |IXB 



below. First, however, let us argue qualitatively where the difference between parafermions (qubits) and fermions, 
bosons originates from. An example will illuminate this. For the case of bosons and fermions, [b\b2,b\b3] ~ b\b3 
and [/1/2, 72/3] = fifs- But for parafermions, [a|a2,a2a3] = a|a3(l — 2n2). (An easy way to check this, without 
explicitly calculating the commutator, is to use the mapping to fermions: fjfi+i ^ ajoi+i and the Jordan- Wigner 
transformation /j — > aiSi.) Thus the difference is that bosons and fermions preserve locality, but parafermions do not. 

Similarly, we can consider additional bilinear operators. For fermions, if we also have fifj and fjfj, the group 
is SO{2N), which is too small to be interesting for QC. In fact this is a model of non-interacting fermions: there 
exists a canonical transformation to a sum of quadratic terms each of which acts only on a single mode (see also 
Refs. |6|, ^ , ^0|j3^j37[| ) . For bosons, if wc include bibj and bjbl, the group generated is the N{2N + l)-parameter 
symplectic group Sp{2N, R) which is non-compact, implying that it has no finite dimensional irreps ||l9[| . If we 
further include the set of annihilation and creation operators bi,bl together with the identity operator /, the set 
{I,bi,bl,bibj,b^jbl,bjbi} generates the semidirect-product group N{N) (g) Sp{2N,R), where N{N) is the Heisenberg 
group, with {N + 1){2N + 1) generators (Ref. Ch.20]). This is therefore still too small to be interesting for 
universal QC. In fact, this is exactly the reason why linear optics by itself is insufficient for universal QC. The 
situation does not change even after introduction of the displacement operators Di{a) — exp(a6|— a*bi) [p^ , since 
Di{a) £ N{N)^ Sp{2N,R). 

The way to universality [with respect to SU{2^)] is to introduce nonlinear operations such as a Kerr nonlinearity 
p8| , self- interaction or conditional measurements A Kerr nonlinearity is a two-qubit interaction of the 

form n^n'j (where i and j are different modes), which directly provides a CPHASE gate. To see this, consider a 

dual-rail encoding [^: Suppose that one qubit is encoded into |0) — b\\0), |1) — b\\0), while a second qubit is 

l\o),\i) ^bl\ 

|00) = blb\ |0) , |01) = blbl |0> 
\10) =blb\\0), |ll)=5l4|0> 



encoded into |0) = 63 |0), |1) = 64 |0) (|0) is the vacuum state). The two-qubit states are 



(This is related to Case 3 of section VI, where we showed how to make qubits from bosons.) It is then simple to 
verify that cyiYi{—i'!m\n\) acts exactly as a CPHASE gate, i.e., it is represented by the matrix diag(l, 1, 1, —1) in this 
two-qubit basis. Here we wish to point out that a recently introduced alternative to a Kerr nonlinearity | |3l[ |, namely 
the self-interaction {n\Y , is in fact closely related to the Kerr nonlinearity. Thus methods developed to use one of 
these non-linear interactions can be transported to the other. Let us demonstrate this point by giving a simple circuit 
to show how one interaction simulates the other. We start with the operator identity 

exp(0(at6 — b^a))b^ exp(— 0(a^5 — b^ a)) — cos 06^ + sin^a^, 

which can be proved directly from the Baker-Hausdorff formula 

e-^^Be^^ = B- a[A, B] + ^ [A, [A, B]] + ^[A, [A, [A, B]]] + ... (12) 
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Using the latter identity it is then simple to verify the following identity, which holds on the two-quhit subspace above, 

exp(-i7rn°n^) = &L^{--{b\W - 6I62)) exp(-z7r-^^^ -)&^\,{-{h\h^ - h\h-2)). 

This is an exact 3-gate simulation of the Kerr CPHASE gate in terms of the self-interaction. The simulation uses 
the linear bosonic operators and the local energies n,^ in order to unitarily rotate the self-interaction terms 
(n|)^ + ("-4)^ to a Kerr interaction. 



VIII. FLUCTUATIONS IN PARAFERMION NUMBER AT FINITE TEMPERATURE 



So far we have not really made use of the full power of the Fock space representation, which allows to consider 
the case of fluctuating particle number. The quantum statistics of parafermions is determined by their commutation 
relations, like fermions (Fermi-Dirac statistics) and bosons (Bose-Einstein statistics). A simple case to consider is 
that of non-interacting parafermions. The Fermi-Dirac distribution for an ideal Fermi gas is derived using only the 
restriction that no more than a single fermion can occupy a given mode ]39[ | . Hence the statistics of non- interacting 
parafermions is clearly the same as that of non- interacting fermions. 

Fluctuations in particle number will be a result of interaction of the system with an external bath, which imposes 
a chemical potential /i (essentially the gradient of the particle flow). As a simple example, consider the following 
system-bath interaction Hamiltonian: 

N N 

^^(2n,-l)®i3^ (13) 

4=1 1=1 

where are bath operators. To further simplify things assume the bath is treated classically, i.e., are positive c- 
numbers. With this Hamiltonian, one can study the fluctuations of parafermions under finite temperature T. Mapping 
from the well-know result for a non-interacting Fermi gas it then follows that the average occupation for the i^^ 
qubit site is 

('^^) = g(2i3^--M)/feT_^^ 

where k is Boltzman's constant. This is the average value of the qubit-"spin" (whether it is |0) or |1)). Keeping the 
chemical potential fi fixed, in the limit of T ^ we find that (n;) ^ 1 if < /i, but (rii) — > if > ^. Thus, as 
expected, it is essential to keep the interaction with the bath weak (compared to /i) to prevent fluctuations in qubit 
"orientation" at low temperatures. At finite T we find (n^) < 1, meaning that some fluctuation is unpreventable. Of 
course, our model is very naive, and the picture is modified when qubit interactions are taken into account. However, 
it should be clear that a Fock space description of qubits, i.e., in terms of parafermions, could be valuable in studying 
qubit statistics at finite temperatures. 



IX. UNIVERSALITY OF EXCHANGE- TYPE HAMILTONIANS 



In this final section we conclude with an application of the formalism we developed above to the study of the 
universality power of Hamiltonians. We have considered this question in detail before for general exchange-type 
Hamiltonians (isotropic and anisotropic) jl^jl^. We first briefly review the universality classification of various 
physically relevant bilinear Hamiltonians. It will be seen that while in certain cases the Hamiltonian is not sufficiently 
powerful to be universal with respect to U{2^), it is universal with respect to a subgroup. As mentioned in Section 



III 



this result requires the use of encoding of physical qubits into logical qubits pl[-p4| . We then consider in detail 
the representative example of the XY model, where we give a new result about universality (in fact, the lack thereof) 
in the case of nearest-neighbor-only interactions. We then present new results about codes with higher rates than 
considered in Refs. ||l^,|I^. For simplicity we revert when convenient to the Pauli matrix notation in this section, 
which is more familiar to practitioners of QC. 
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A. Classification of Bilinear Hamiltonians 



The most general bilinear Hamiltonian for a qubit system is 

H{t) ^Ho + V + F^J2 + E + ^' (14) 

i i<j 

where Hq is the qubit energy term, the interaction between qubits i and j is: 

a.l3—x.y,z 

and the external single-qubit operations are: 

i 

Recall the "standard" result about universal quantum computation: The group U{2^) on N qubits can be generated 
using arbitrary single qubit gates and a non-trivial two-qubit entangling gate such as CNOT p^ . The general 
Hamiltonian H{t) can generate such a universal gate set, e.g., as follows: Suppose there are controllable (t| and erf 
terms. Then af can be generated using Euler angles: 

af = expi-i-KCT- /4:)a^ exp{i'Ka^ /4:). 

This is an instance of a simple but extremely useful result: let A and B be anticommuting hermitian operators where 
= I {I is the identity matrix). Then, using Ue^W = e^^^ (U is unitary, V is arbitrary): 

o exp{i6B) = exp{—iAip) exp{i9B) exp{iALp) 
exp{-i9B) if = 7r/2 



exp[id{iAB)] if = 7r/4 ' (^^^ 

One can also derive these relations for su(2) angular momentum operators, without assuming that {A, B} — and 
A^ = /. Let Jx and be generators of su{2). Then, using the Baker-Hausdorff relation Eq. (p^, and [J^, J^] — iJy- 

exp{—iipJz)Jx exp{i(pJz) = Jx cos 1^9 + Jy svaip. 

From here follows, using C/e^C/^ — e^^^^ again: 

C'j^ o exp{i9Jx) = exp{i0{Jx cos ip + Jy sin 93)), 

and Eq. (|l|) can be verified, with ip —^ 2ip. 

Different QC proposals usually have different two-qubit interactions. Typical types include erf cr,f_|_]^, cr*'(T^_^j^ (or 
crfaf_^i), afaf_^_i + crfcrf(.i (XY model), and ifi ■ ifj (Heisenberg model). It is simple to show that they can all be 
transformed into a common canonical form cr|^crf ij, using a few unitary transformation. The term afa^_^_l can be used 
to generate CPHASE and from there, CNOT E.g., the XY term can first be transformed into crfafj^i using 

Euler angle rotations about erf, which flips the sign of the cr^af_^^ term: 



exp 



Y(«+i+'Tf'^r+i) 



L^L o exp 



2 

= exp(i6i(Tf crf+i), 

which can subsequently be transformed into the canonical form using another Euler angle rotation: 

C^y^^y^^ o cr- (T^j^-^ = cr^ cr,-+i, 

where using ^^0^,-,] = we have abbreviated C'"i^ o C'^i^ as C'^/'t . The method of Euler angle rotations as 

O L 2 ' Z + iJ + 1 

applied here is also known as "selective recoupling" in the NMR literature . 

Not all QC proposals have an interaction Hamiltonian that appears to be of the form Vij. E.g., the ion-trap proposal 
[^2| looks quite different since it involves interactions between ions mediated by a phonon. The interaction between 
the i*'' ion and the phonon has the form (t,^6^ + crfh. This is nevertheless equivalent to an XY model, since: 



i+l 



Therefore in many cases it suffices to study the interaction a^afJ^l. 
Let us now consider a number of more restricted models: 
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1. No external single-qubit operations 
If F = then the nearest- neighbor set {erf , crfcrf^^ , erf crf^^, o-f^_^CT| } is stih universal, since 

This is the case when Hq is controUable. More physically, the set {erf , Ir^i • Vi+i, (Vi x = (JiCr^j^T^ — crf^^^crf } 

is also universal, where 'a' = {a^ , a^). The term If'i x "t^i+i is an anisotropic (Dzyaloshinskii-Moriya) interaction 
which arises, e.g., in quantum dots in the presence of spin-orbit coupling |p7| , ^ -^ . 



2. No external single-qubit operations and Ho uncontrollable 

If F = and Hq is not controllable, then the nearest-neighbor set {<Jiaf_^_i,afaf_^_i,afaf^^,afaf^-i^,a^a^_^_-^] is 
universal, meaning that the interaction term V by itself is universal. One way to see this is to map the set to 
parafermionic operators and note that it overlaps with the set that generates the parafermionic C/(2^) (Table 



3. Scalar anisotropic exchange-type interactions 

Consider the case J^"'^ — J"jSai3 (denoting V hyV), which amounts to limiting the Hamiltonian to scalar anisotropic 
exchange- type interactions. Using Eq. (0) we then arrive at the second-quantized form 



i 



F' = ^ Ay (ojOj -1- a|a]) -f- Jij{a\aj + a]ai) -f AJ^^UiUj (16) 

where 



^,=£,;+ ( I , ^ (/f - ), 



and we dropped a constant energy term. 

V is the so-called XYZ model of solid-state physics. Considering the structure of V and the classification of 



operator algebras we carried out in Sections [II, IV, it should be clear that some immediate conclusions can be drawn 
about the universality power of this Hamiltonian. The full Hamiltonian + V"' + F contains the generators of the 
parafermionic t/(2^) (Table ^, so it is universal. On the other hand, without external single qubit operations F = 0, 
whence [i?o + — 0, so -V V' sSAp, i.e., preserves parity. This immediately implies that the XYZ model (even 
with Hq) is by itself not universal. However, it can be made universal by encoding logical qubits into several (two are 
in fact sufficient) physical qubits [|l^. The elimination of single qubit operations {F = 0) can be quite useful, since 
typically single and two-qubit operations involve very different constraints. In some cases single-qubit operations can 
be very difficult to implement (sec [|l^,|l8|,^ and references therein for extensive discussions of this point). 



B. XY Model 



Consider now the XY model, which is defined by 

VxY ='^Jrj (4 + a] ) . (17) 

i<j 
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It is relevant to a number of proposals for quantum computing, inclu ding quantum Hall systems E3,E6i, quantum 



dots in microcavities quantum dots coupled by exciton exchange |47], and atoms in microcavities|l48| . Let us 



summarize what is currently known about quantum computational universality of this model. 

• In Ref. it was shown that the XY model with nearest neighbor interactions only, together with single qubit 
operations, is universal. 

• In Ref. it was argued that the XY model is universal without single qubit operations, provided these gates 
can be applied between nearest-neighbor and next-nearest-neighbor pairs of qubits. This involved encoding 
a logical qutrit into three physical qubits: |0l) = |001), = |010), \2l) — |100). We reconsider this in 
subsection [X D| below in the context of the XXZ model (but using the methods of the results arc valid also 



for the XY model). 

In Ref. |l|] we showed that the XY model is universal using only nearest and next-nearest neighbor {Ji^i^2) 
interactions, together with single qubit cr^ terms. This too involved an encoding, of a logical qubit into two 
physical qubits: |0l) = |01), |1l) = |10). Two comments are in order about this result: first, next-nearest 



neighbor interactions can be nearest neighbor in 2D (e.g., in an hexagonal array); second, unlike |32|, we did not 



assume the CTz terms to be controllable, i.e., there is no individual control over Si [Eq. (|14|)]. A similar model is 



treated in subsection IX C below 



The question now arises: Is the XY model universal with nearest-neighbor interactions only? We prove that it is 
not. 

The nearest-neighbor XY model in its parafermionic form is 



JV N 



Consulting Table |l|, we see that H can only generate SU{N), which is clearly too small even for encoded quantum 
computation. 



C. Antisymmetric XY Model 



To illustrate the idea of encoding for universality, let us briefly consider the "antisymmetric XY model" : 

V.XY = E J^atay + Jffafa^. (18) 

Here J^J^ and JfJ' are real. We encode a logical qubit into pairs of nearest-neighbor physical qubits. Letting 

= Jif ~ ^ = Jij' + Jij'^ = £2m-l " £2™, (19) 

using the compact notation = ■nm-\.2m-, and assuming that interactions are on only inside pairs of qubits encoding 
one qubit, we find for the Hamiltonian H = Hq + V^xy'- 

N/2 

H.XY = Yl {J^^'n + 4^Kn) + {^mTl + r„^,) , (20) 
in—l 

where the T and R operators were defined in Eqs. (|),(0). Since the T and R operators form commuting sl{2) algebras, 
the Hilbert space splits into two independent computational subspaces. The R operators conserve parity, so that an 
appropriate encoding in the axially symmetric case (A^ = 0), using standard qubit notation, is |0l) = |00) and 
|li) = |11). On the other hand, the T operators preserve particle number, so that if Jm = (axially antisymmetric 
case) the encoding is |0l) = |01), |1l) = |10). In both cases control over the pair of parameters {Jm,e+} (or 
{l^rm^mY) Sufficient for the implementation of the single-encoded-qubit SUm(2) group (the subscript m refers to 
the m^^ logical/encoded qubit). 

Logic operations between encoded qubits require the "encoded selective recoupling" method introduced in Ref. . 
Consider the "axially antisymmetric qubit" |0l) — |01), |1l) = |10). First, note that using Eq. (p^: 
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C^Tfl ° ^23 = i(^l(^2Ti3- (21) 

Now assume we can control A13; then: 

Cr'i ° {Ct^, o TI3) = al {al ~a{)l2. (22) 
Since a\a2 is constant on the code subspace it can be ignored. On the other hand, erf erf acts as — r|: 

|0l)i|0l)2 = 101)12101)34 -101)12101)34 - -|0l)i|0l)2, (23) 

and similarly for the other three combinations: |Oi)|U) -> |0l)|U), |1l)|0l) ^ |1l)|0l), |1l)|1l) ^ -|1l)|1l)- I.e., 
(T|cr| acts as an encoded cr^ ® . This establishes universal encoded computation in the antisymmetric XY model. 



D. Codes with Higher Rates 

The encoding of one logical qubit into two physical qubits is not very efficient. Can we do better? I.e., can we 
perform encoded universal QC on codes with a rate (no. of logical qubit to no. of physical qubits) that is greater 
than 1/2? We will show how in the case of the XXZ model, defined as H = Ho + Hxxz, where 



When surface and interface effects are taken into account, the XY-examples of QC proposals [ ^2|j4^ -j4^, as well as 
the Heisenberg examples ||l^,^,Q , are better described by the axially symmetric XXZ model. Additional sources of 
non-zero J^j in the XY-examples can be second-order effects (e.g., virtual cavity-photon generation without spin-flips 
p2| ) . A natural XXZ-example is that of electrons on helium |3lj3^ . 

First, note that the code used in the XY model, |0l) =^01)^ |1l) = |10), is applicable here as well: Tf^ = 
^{o'fcTj +cr^(Tj) preserves particle-number, and serves as an encoded a^; erf terms from Hq serve as encoded cr^, and 
ffcrf^j apphed to physical qubits belonging to different encoded qubits acts as encoded i^) . 

In the general encoding case we consider a block of N qubits where codewords are computational basis states 
(bitstrings of O's and I's): {q^iNa) \0)}a, where a — {ai} and ai can be or 1, while Na = 0...N. A code-subspace 
C(-/V, n) will be defined by having a fixed number n of I's (i.e., of parafermions). Thus there are 



dx.n = dim[C(iV, n)] 



codewords in a subspace. Examples are considered below. Note that these subspaces are decoherence-free under the 
process of collective dephasing p^ , and have been analyzed extensively in this context in Ref. (52j. Figure 1 in 
Ref. provides a nice graphical illustration of the C{N,n) subspaces. Since the decoherence-avoidance properties 
of the codes we consider here have been extensively discussed befor e l5l| , |5^ , and even implemented experimentally 
we do not address this issue here. We further note that Ref. [|S2| provided an in-principle proof that universal 
encoded QC is possible on all subspaces C{N,n) independently. However, this proof had several shortcomings: (i) it 
used a short-time approximation, (ii) it did not make explicit contact with physically realizable Hamiltonians, (iii), it 
proceeded by induction, and thus did not explicitly provide an efficient algorithm for universal QC. We remedy all 
these shortcomings here, i.e., we (i) use only finite-time operations, (ii) use only the XXZ Hamiltonian, (iii) provide 
an efficient algorithm that scales polynomially in N. 

We need a measure that captures how efficient a C{N,n) code is. If there are d codewords, supported over N 
p-dimensional objects (p = 2 is the case of bits), and information is measured in units of q, then we define the rate of 
the code as 

log d 

r{d,p,q) - 



logg ' 

The traditional definition for qubits is recovered by setting p = q = 2, i.e., the rate of a code is the ratio of the number 
of logical qubits log2 d to the number of physical qubits A^, which in our case becomes: 

r = S{e) (24) 
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where e = k/N , 

5(e)--elog2e-(l-e)log2(l-e) 

is the Shannon entropy, and we have used the Sthding formula logx! « xlogx — x. Since S{\) — 1 the code has 
a rate that is asymptotically unity for the "symmetric subspace" C{N,N/2), where the number of I's equals the 
number of O's in each computational basis state. However, we will not in fact attempt to encode log2 dAr,rt logical 
qubits in the subspace C{N,n), since the subspace does not have a natural tensor product structure. Instead we will 
consider C{N,n) as a subspace encoding a qudit, where d — dN,n- Using the generalized definition of a rate above, 
and measuring information in units of d so that each subspace encodes one unit of information, the rate of such a code 
is r = 1^^^. This, however, exactly coincides with r of Eq. (|2^). Therefore we see that the advantage of working 
with the symmetric subspace C{N,N/2) in the limit of large N is that its rate approaches unity. 

Before embarking on the general analysis, let us note that for an encoding of one logical qubit into TV physical qubits, 
there is a simple construction in terms of parafermionic operators: Qa,f3{N), pi^), and [Q]^ /jC-^); Qa,/3(-^)] (which 
is a function of parafermion number), form an su{2) algebra in the basis \0l) = qti-^a) |0) and — qjj{N ~ Na) |0). 
E.g., for iV = 2 there are two cases: the sets {0102, a,2a\, fii + h2 — 1} and {a\a2, 0201, fii — 71,2}, with corresponding 
bases \0l) = |0), = a\al |0) and \0l) = a\ |0), — al \^)- These two encodings are universal (in the sense of 
blocks of N physical qubits) when only Hq and V' are controllable [Eq. ([l6|)]. 

Let us now move on to the general subspace case, starting with an example. 



1. Encoded Operations: Example 

Consider C(3,l) — Span{|0) = |001),|1) = |010),|2) = |100)}, i.e., an encoding of a logical qutrit into 3 physical 
qubits, as in Ref. |^^. Let us count qubits as i = 0, iV — 1. Our first task is to show how to generate su(3) on this 
subspace. It is simple to check that TJ^lOOl) = 0, To^ilOlO) = |100), T^illOO) = |010), and in total 




To^i = 1 = |1)(2| + |2)(1| ^Xi2 

\o 1 0; 

where the notation X12 denotes a operation between states |1) = |010) and |2) = |100). Similarly it is simple 
to check that = ^01 and Tq2 = Xo2- Further, using T^^ = i (erf — ""l) we have: Tq\ = Z12, Tfj = ^01, and 
Tq2 = Z02, where Z12 denotes a operation between states |1) and |2), etc. Therefore each pair {Tj^, T?} generates 
an encoded su{2). But in the sense of generating, su{N) is a sum of overlapping su(2)'s so using just the nearest 
neighbor interactions {Tqi,T^i,Ti2,Ti2} we can generate all of sm(3) on C(3, 1). Note that [Xoi,Xi2] — ilo2, so that 
su{2) between states |0), |2) can in fact be generated using Tf^'s alone, without Tij's. This conclusion clearly holds 
for the generation of all of su{3) on C(3, 1), as first pointed out in Ref. p3[ . 

Next, we need to show how to implement encoded logical operations between two C(3, 1) code subspaces. Let us 
number the qubits as i = 0, 1, 2 for the first block, i = 3,4,5 for the second block. Consider the effect of turning on 
JI3, i.e., consider the action of tT|tT| on the tensor product space C(3, 1) ® C(3, 1). The operator cr|cr| is represented 
by a 9 -dimensional diagonal matrix on this space, which is easily found to have the following form in the ordered 
basis {|0)(^ |0),|0)O |1),...,|2)(^ |2)}: 

a^al = diag(-l, 1, 1, -1, 1, 1, 1,-1,-1) = diag(-l, 1, 1) ® diag(l, 1, -1). 

E.g., cr|cr||2) |2) = crfcrfllOO) (g) |100) = (+|100)) «) (-|100)) = -|2) (g) |2), which explains the -1 in the 9"^ position 
in the diagonal matrix above. The important point is that (T2(t^ acts as a tensor product operator on C(3, 1) (g)C(3, 1), 
which puts a relative phase between the basis states of each C(3, 1) factor. This means that cr|cr| acts as an "su(3) 
-like" a'^ (g) on C(3, 1) (g) C(3, 1). (It is an "sw(3)-like" tr'^ (g since for su{2) a"" = diag(l, -1) and here we have 
instead diag(— 1,1,1) and diag(l, 1, — 1).) It is well known that the CPHASE gate can be generated from the 
Hamiltonian (g . The same holds here, so that we can generate a CPHASE gate between two C(3, 1) subspaces 
by simply turning on a nearest-neighbor interaction between the last qubit in the first block and the first qubit in the 
second block. 

With this example in mind we can move on to the general case. 
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2. Encoded Operations: General Subspace Case 

Let us now consider the case of a general subspace C{N,n). We can enumerate the codewords as {|0), \dN.n)} 
where |0) = |0, 01, 1) etc., to |c?7v,n) — |1, 10, 0), where there are TV qubits in total and n I's in each 
codeword. Consider a fixed nearest-neighbor pair of qubits at positions i,z + 1, and the action of Tfj_|_]^, r/j_|_]^. The 
four possibilities for qubit values at these positions are {00,01, 10, 11}. Now consider a pair of codewords such 
that \t) has 01 in the + l positions while \t') has 10 in the + l positions, and they are identical everywhere else. 
We can always find such a pair by definition of C{N,n). The action of Tf T/j_,_j on is to generate su[2) 

between them, just as shown in the case of C(3, 1) above. On the other hand the action of T^^j^^, T.^^j^^ in the case of 
00 or 11 in the z,z + 1 positions is to annihilate all corresponding codewords (which are anyhow outside of the given 
C(iV, n) subspace). This null action means that, when exponentiated, Tf^_^_^, T^^_^^ act as identity on these codewords. 
Therefore the action of ?"fi+i, Ti^i+i is precisely to generate su{2) between \t),\t'), and nothing more. Denote this by 

su{2)^l'}j^i. Let us now keep the 01 and 10 at positions i,i + 1 fixed, and vary all other N — 2 positions in \t),\t'), 

subject to the constraint of n 1 's, and in the same manner in both We then run over K ~ (^rf) codewords, 

and Ti^i+i generate su{2) between each pair of new \t),\t'). Denote these by sm(2) k = 1..K. By further 

letting 1 = 0, ..N — 2 we generate N —1 overlapping su(2)'s. These sm(2)'s can be connected by swaps so that we can 
generate all su(2)^^\ k — 1..K, i < j. We thus have a total of (^Zi) (^) su(2)'s. To generate the entire su{dN_n) 
we need no more than d^^n = (^) overlapping su(2)'s. Since (^Zf) (^) / (^) = ^n{N — n)>l, we have more than 
enough overlapping sm(2)'s, and su{dN^n) can be generated. 

What is left is to show that we can perform a controlled operation between two C{N,n) subspaces. To do so we 
again use the nearest-neighbor interaction (t|^_]^(t^, where the first factor {aff_i) acts on the last qubit {N — 1) of the 
first C{N,n) subspace, and the second factor (u^) acts on the first qubit (TV) of the second C{N,n) subspace. Now 
let us sort the codewords in the two subspaces in an identical manner, e.g., by increasing binary value. Then consider 
the action of ct^.^ct^ on the resulting ordered basis {|0) (g) |0), |0) (g) |1), |dAr,„) (S) \dN,n)}- This action generates a 
representation of cr^_]^(T^ by a dj^^n x dj^^n diagonal matrix. As in the C(3, 1) case considered above, this matrix is 
actually a tensor product of an "su((iAr_„)-like" cr^ (g) on C{N, n) ® C{N, n). It is simple to determine the form of 
these two (different) cr^'s. For the codewords belonging to the left C{N, n) factor write down a -1-1 (—1) for each (1) 
in the 7V'^ position. These numbers are the diagonal entries of the left "su(dAr.„)-like" cr^ factor. Similarly, for the 
codewords belonging to the right C{N,n) factor write down a -1-1 (—1) for each (1) in the A'^ -I- l"^ position. These 
numbers are the diagonal entries of the right "su((iAr,„)-like" factor. Since each such "su((iAr,„)-like" puts relative 
phases between the basis states of C{N,n), the action of ct^^^^ct^ is that of a generalized CPHASE between the two 
code subspaces. This is sufficient together with su(dN „) on each block to perform universal quantum computation 
|6|. 

X. CONCLUSIONS 

The standard quantum information-theoretic approach to qubits and operations on qubits, emphasizes qubits as 
vectors in a Hilbert space and operations as transformations of these vectors Q . This is the point of view of the first- 
quantized formulation of quantum mechanics. An alternative, mathematically equivalent, point of view is the Fock 
space, second- quantized formulation of quantum mechanics, which emphasizes the particle-like nature of quantum 
states. Qubit up/down states are replaced by qubit presence/ absence, while rotations are replaced by operators 
that count or change particle occupation numbers. The mapping of qubits to parafermions discussed in this paper 
is a mapping between these first and second quantized formulations. It proved to be a useful tool in studying the 
connection between qubits, bosons and fermions, in analyzing the algebraic structure of qubit Hamiltonians, and in 
studying related quantum computational universality questions. In particular, it allowed us to classify subalgebras of 
fermion, boson, and qubit operators and decide their power for quantum computational universality. These results 
are relevant for physical implementation of quantum computers: a physical N -qubit system comes equipped with 
a given Hamiltonian, which generates a subalgebra of su{2^). It is important to know whether this Hamiltonian is 
by itself universal or needs to be supplemented with additional operations, or whether one needs to encode physical 
qubits into logical qubits in order to attain universality. Our classification settles this question for many subalgebras 
of physical interest. 

Another potential advantage of the parafermionic approach, as a second-quantized formalism for qubits, lies in its 
ability to naturally deal with a "qubit-field" , i.e., situations where the qubit number is not a conserved quantity. This 
is certainly a concern for optical and various solid-state quantum computer implementations. We leave the study of 
a qubit field theory as an open area for future explorations. 
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group 


fermions 


parafcrmions 


U{2") 






SO{2N + 1) 


/l/.,/«/.,/.,h.c. 


Qjj^ SiSj dj ^ diSi Sj dj , diSi , ll.C. 


SO{2N) 


h-c. 


Si S j dj ^ diSi S j dj , h .c. 


U{N) 


flfj 


Si Sj dj 



TABLE I. Infinitesimal generators (h.c.=hermitian conjugate). 



group 


formions 


parafermions 


17(2") 


/»S/,///.+i,h.c. 


Oi, ajai+i, h.c. 


SO{2N + 1) 


/i,h.c. 


OiSi, h.c. 


SO{2N) 


//./'i+l,,f»,fi+l,h.C. 


alai+i,aiai+i, h.c. 


SU{N) 


///»+i,h.c. 


ajfli+i, h.c. 



TABLE II. Hamiltonian generators. 
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